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Comparison theorems are established for the Dirac and Klein-Gordon equations. We suppose that 
V'^^\r) and V^'^\r) are two real attractive central potentials in d dimensions that support discrete 
Dirac eigenvalues S^^'^ and E^^^^ ■ We prove that if 1/'^' (r) < V^'^' (r) , then each of the corresponding 

discrete eigenvalue pairs is ordered e\^^^ < E^^\. This result generalizes an earlier more restrictive 
theorem that required the wave functions to^ae node free. For the the Klein-Gordon equation, 
\ similar reasoning also leads to a comparison theorem provided in this case that the potentials are 

^SJ ■ negative and the eigenvalues are positive. 
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^ ! I. INTRODUCTION 

■ The comparison theorem of quantum mechanics states that if two real potentials are ordered, V^^^{r) < V^-^^ (r) , then 
I .... (1) (2) 

. each corresponding pair of eigenvalues is ordered E^J < E^J . For non- relativistic problems whose Hamiltonians are 

■ bounded below, this theorem is an immediate consequence of the variational characterization of the discrete spectrum. 
I— I, By contrast, since the energy operators for relativistic problems are not bounded below, a simple comparison theorem 

is unexpected since the usual variational arguments are not applicable l'-3]. However, by using reasoning developed 
originally for Schrodinger problems in which the graphs of the comparison potentials cross each other in a controlled 
way and consequently variational inequalities are no longer available, it was possible to derive a limited comparison 
theorem Q valid for node- free Dirac eigenstates; this was subsequently extended to Dirac problems in d dimension 
In a more recent development 7] , it was shown for the Dirac problem that if an attractive potential V{r, a) is monotone 
in a potential parameter a, then each corresponding Dirac eigenvalue E{a) is monotone in this parameter. Both types 
of result have also been proved for potentials in the Klein-Gordon equation 0], with the added restriction in this 

■ case that the eigenvalues considered are positive. We shall see in this paper that the monotonic spectral behavior 
' induced by attractive central potentials V{r, a) that are monotonic in a parameter a leads to a general comparison 

theorem for the Dirac problem, and to a comparison theorem restricted to negative potentials and positive eigenvalues 
for the Klein-Gordon equation. These results are established below in sections H and HI respectively. In section IV 
we consider an example of comparison approximations to the Dirac spectrum of a screened-Coulomb potential: in 
' particular this establishes a general energy-bound conjecture posed 25 years ago in Ref. [91] and later proved to be 
true for node-free states in Ref. Q- 

The potential V{r) discussed in this paper is sometimes called a 'vector potential' since it enters the problem as the 
time component of a 4-vector: it is this potential that for Hydrogen-like atomic models is taken to be V{r) — —Za/r. 
A 'scalar potential' S{r) may also be considered, as a variable term added to the mass [l0|. However, throughout 
the present paper, for both Dirac and Klein-Gordon problems, we shall assume that S = and that the mass m is 
constant. 
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II. THE DIRAC EQUATION 



For a central potential V(r) in d dimensions the Dirac equation can be written in natural units ?i = c = 1 as 

i—^H^, where H = 2_^asPs + mP + V, (1) 
ut 

s—l 

where m is the mass of the particle, and {a^} and /? are Dirac matrices, which satisfy anti-commutation relations; 
the identity matrix is implied after the potential V . For stationary states, algebraic calculations in a suitable basis 
[m lead to a pair of first-order linear differential equations in two radial functions {^i(r), -02 ('')}, where r = |r|. For 
d > 1, these functions vanish at r = 0, and, for bound states, they may be normalized by the relation 

oo 

(V'i,Vi) + (^2>2) = J{i>l{r)+ijl{r))dr = l. (2) 



We use inner products without the radial measure factor ' because the factor r 2 is already built in to each 
radial function. Thus the radial functions vanish at r = and satisfy the coupled equations 

E4,i = (V + m)i;i + {-d + kd/r)^2 (3) 
Eij2 = {d + kd/r)^i + {V-m)i^2, (4) 

where fci = 0, fc^ = T{j + ^^), d > 1, t ^ ±1, and E = Ek^v is the eigenvalue corresponding to the state with 
J' = 0, 1, 2, . . . nodes in the upper radial function ipiir). We note that the variable t is sometimes written u, as, for 
example in the book by Messiah [l2|, and the radial functions are often written ipi — G and -02 = F, as in the book 
by Greiner We shall assume that for each attractive central potential l/considered, there is a discrete eigenvalue 
E = Ek^y and that Eqs.(I31l3]) are the eigenequations for the corresponding radial eigenstates. In this paper we shall 
present the problem explicitly for the cases d> \. Similar arguments go through for the case d= 1: in this case ki —Q, 
the states can be classified as even or odd, and the normalization ([2]) becomes instead [tpKx) + 01 (^)) dx — 1. 
We first suppose that the potential V — V{r, a) depends smoothly on a parameter a, and we re-state a theorem from 
Ref. H: 

Theorem 1 The real attractive central potential V{r,a) depends smoothly on the parameter a, and E[a) ~ Ek^uio-) 
is a corresponding discrete Dirac eigenvalue. Then: 

dV/da > E'{a) > and dV/da < ^ E\a) < 0. (5) 

The principle result of this section of the present paper is the proof of the following general comparison theorem. 

Theorem 2 Suppose that E^}^^^ and E^^J^ are Dirac eigenvalues corresponding to two distinct attractive central 
potentials V'^^''{r) and V'^^''{r). Then: 

<y (2) (,) ^ 4;,i<4'l. (6) 



Proof of Theorem 2: The theorem may be proved by keeping the given real comparison potentials V'^^\r) and 
V^^\r) fixed and defining the one-parameter family of potentials V{r,a) by 

Vir,a) = V'^^'>{r)+a(v'^^\r)-V^^\r)y aG[0,l]. (7) 

We now consider one of the discrete eigenvalues, Ek^^{a) = E{a), of the Dirac equation with potential V{r,a). In 
terms of this energy function, the comparison eigenvalues of the theorem are given by eI^^^ = E{0) and E^^^^ = E{1). 
Clearly, since V^'^Hr) - V^^^\r) > 0, we have dV{r,a)/da = V^'^\r) - V^^^^r) > 0. Thus, by Theorem 1 we have 
E'{a) > 0. The fact that E{a) is increasing in turn implies that E{0) < E{1), which result establishes the theorem. 



□ 



3 



III. THE KLEIN-GORDON EQUATION 



The Klein-Gordon equation for an attractive central potential V{r, a) depending on a parameter a is given by 

(-A + to2)V' = (i; - F) V- (8) 

The generalization to arbitrary dimension d > 1 is conveniently allowed for Q by replacing the orbital-angular 
momentum quantum number £ of the three-dimensional problem by Id-, where 

id^e + {d~3)/2. 

The eigenvalues of the problem in d dimensions are labelled by £d and the number — 0,1,2, ... oi nodes in the radial 
eigenfunction '^(r). The radial equation becomes from ([5]) 

- + = ((S,,. - V{r)f - m') ^{r), (9) 

where 

Q = i(2^ + d - l){2£ + d-3)^ id{£d + 1), £^0,1,2,..., d^2,3,... (10) 
In Ref. Q the following theorem was proved: 

Theorem 3 The real attractive central potential V{r, a) < depends smoothly on a parameter a, and E{a) = 
E£^i,{a) > is a corresponding discrete Klein-Gordon eigenvalue. Then: 

dV/da > ^ E'{a) > 0, and dV/da < ^ E'{a) < 0. (11) 



By similar reasoning to that used for the Dirac case, we shall now prove the Klein-Gordon comparison theorem: 

Theorem 4 Suppose that E^^^J^ and Ef'^^ are Klein-Gordon eigenvalues corresponding to two distinct attractive central 
potentials V^^^r) andV^^^r). Then: 

V'^'Hr) <V^^\r) <0 and eIII>0 E^ < e'^H (12) 



Proof of Theorem 4: We suppose that y(i)(r) < F(2)(r) < and we define a one-parameter family of negative 
potentials V{r,a) by 

V{r,a) = V'^^Hr)+a(v'^^\r)-V^^\r)^ , a e [0,1]. (13) 

Clearly, V{r,a) < 0. Meanwhile, dV{r,a)/da = V^^Hr) - V^^^^r) > 0. We also suppose that E^H = E{0) is positive. 
By Theorem 3 we deduce that E'{a) > for a G [0, 1]. Thus we conclude that 

E^ll^Eil)>EiO)^E^,ll, 

which inequality establishes the theorem. 

□ 



IV. ENERGY BOUNDS FOR A SCREENED-COULOMB POTENTIAL 



We revisit an example discussed in Refs. [5|,(9l|, but now we consider the problem in general spatial dimension d > 1, 
and we are able to construct energy bounds valid for every discrete eigenvalue. We first consider the Dirac equation 
for the pure Coulomb problem with potential V{r) = —u/r, where the coupling parameter u — aZ is not too large. 
We write the exact discrete eigenvalues as Dk^v{u) ~ T){u) and they are given [111. [Tsj exactly by 



D{u) = 1 



(14) 
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where 



2 



(15) 



and V = 0,1,2, .. . counts the nodes in the upper radial function 'ipi{r). The principal quantum number n is defined 
in general as 



n = iy+ \kd\ - 



d-3 



(16) 



and the spectroscopic designation 



is provided by the formula 



{s,p,d,...} ^£^{0,1,2,...} 



(17) 

With these conventions, the eigenvalue formula (I14p becomes formally similar to the expression for the three- 
dimensional case: 



D(u) 



1 + 



(18) 



The comparison theorem established in this paper allows us to use the exact eigenvalues of the Coulomb problem, 
for example, to approximate those of a related problem that is nearby in a different sense from what is normally used 
in perturbation theory. Specifically, we consider the Mehta-Patil screened- Coulomb potential V{r) given [TJ] by 



1-(1-1/Z)- 



Ar 



1 + Ar 

where, for example, the appropriate potential parameters for an atomic model in three dimensions are 

v^aZ and A = 0.98aZi 



(19) 



(20) 



The comparison theorem may now be invoked because we can show that, for each value of the parameter t > 0, V{r) 
is bounded above by a shifted-Coulomb potential y^*)(r) of the form 



V{r)<V^''>{r)^^^ + b{t) 



(21) 



This follows by the following argument from envelope theory [I5l - ll7| . If we write V{r) as a transformation V{r) 
g{h{r)) of the hydrogenic potential h{r) = — l/f, then the transformation function g(h) becomes 



g{h) = V 



h + X{l-l/Z) 1 + 



A 



h-X 



(22) 



It follows immediately that the function g{h) is monotone, g'{h) > 0, and concave, g"{h) < 0. This means that g{h) 
lies above its tangents, all of which are of the form a{t)h{r) + b{t), with a{t) > 0, that is to say, shifted attractive 
Coulomb potentials. The coefficients are given by 



a{t) = g'ihit)), bit) = g{h{t)) - h{t)g' {h{t)). 



(23) 



where r = t \s the point of contact between V'^*'\r) and V{r). For each of these upper tangential potentials, the 
spectrum is given exactly (with the aid of Eq. (jl8p ) by the right-hand side of the equation 



E < D{a{t)) +b{t). 



(24) 



This inequality is a consequence of the potential inequality V{r) < V^*'^{r) and Theorem 2. If we optimize over t > 0, 
and effect the change of variable t ^ u ~ g'{h{t)), then the critical point is found to be the same as that of the energy 
function S{u) given by 
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£{u) = D{u) - uD'{u) + V{~l/D'{u)). 

Thus we obtain the following best upper-bound formula, valid for each discrete Dirac eigenvalue, and expressed in 
terms of the pure hydrogenic spectral function D{u) — Dkj^u{u) of Eq. by 

E < min [D{u) - uD'{u) + V{-1/D'{u))] . (25) 

M>0 

This same formula yields an upper energy bound for any potential V{r) that is a monotone increasing and concave 
function of h(r) = —1/r. For potentials that are convex functions of h{r), the same expression provides a lower energy 
bound. Some numerical values for node- free states in dimension d — 3 are given in Ref. and a more extensive 
table of results is given in Ref. @: at that time, the data in the latter table could only be considered as ad hoc 
approximations, although the inequality Eq. (|25|) was proposed then as a conjecture. Since it follows from Theorem 2 
that Eq. (|25p is valid for all the discrete eigenvalues_and for arbitrary spatial dimension d > 1, we may now state, in 
particular, that the values given in Table 1 of Ref. [9] are a priori all upper bounds. We make weaker general claims 
for the Klein-Gordon equation: for the screened-Coulomb example (where g is concave), D{u) would be set to the 
exact Klein-Gordon Coulomb energy, and the energy bound ()25|) is then known, by Theorem 4, to be valid provided 
that the potential is negative and the eigenvalues considered are positive. 

V. CONCLUSION 

If we examine exact relativistic eigenvalue formulas given, for example, in the book by Greiner [l^], we see that, 
whenever a potential depends monotonically on a parameter, the eigenvalues are monotonic in this parameter in 
the same direction. Since the corresponding non-relativistic problems usually behave this way, we are not surprised. 
However, the question immediately arises as to how general such spectral features are in relativistic quantum mechan- 
ics. In this paper, comparison theorems are proved for Dirac and Klein-Gordon problems in which a single particle 
is bound by a central potential. These results are established without the use of any variational arguments. The 
Klein-Gordon result which we have been able to obtain is limited to negative potentials and positive eigenvalues. In 
the Dirac case, there is no limitation save that the comparison potentials each yield discrete eigenvalues in the same 
angular- momentum sector and with the same number of radial nodes. The illustration, involving shifted- Coulomb 
potential that are upper bounds to a screened-Coulomb potential, shows how the existence of a comparison theo- 
rem immediately leads to spectral approximations. For many relativistic problems, this allows the type of spectral 
reasoning and estimates that are commonly employed in non-relativistic quantum mechanics. 
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